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Abstract — Recently, three Sphere Decoding (SD) algorithms 
were proposed for Spatial Modulation (SM) scheme which 
focus on reducing the transmit-, receive-, and both transmit 
and receive-search spaces at the receiver and were termed as 
Receiver-centric SD (Rx-SD), Transmitter-centric SD (Tx-SD), 
and Combined SD (C-SD) detectors, respectively. The Tx-SD 
detector was proposed for systems with Nt < Nr, where Nt 
and Nr are the number of transmit and receive antennas of the 
system. In this paper, we show that the existing Tx-SD detector is 
not limited to systems with Nt < Nr but can be used with systems 
Nr < Nt < 2Nr 1 as well. We refer to this detector as the 
Extended Tx-SD (E-Tx-SD) detector. Further, we propose an E- 
Tx-SD based detection scheme for SM systems with arbitrary Nt 
by exploiting the Inter-Channel Interference (ICI) free property 
of the SM systems. We show with our simulation results that 
the proposed detectors are ML-optimal and offer significantly 
reduced complexity. 

Index Terms — Sphere decoding, ML decoding, spatial modu- 
lation, space-time shift keying, complexity. 

I. Introduction 

Spatial Modulation (SM) Q, El is a recently devel- 
oped low-complexity Multiple-Input Multiple-Output (MIMO) 
scheme that exploits the channel for information transmission 
in an unprecedented fashion. Specifically, the information 
bitstream is divided into blocks of length log2(A^tA/) bits, and 
in each block, logj (M) bits select a symbol s from A/-ary 
signal set (such as A/-QAM or -PSK), and logj (Nt) bits select 
an antenna out of Nt transmit antennas for the transmission 
of the symbol s. The throughput achieved by this scheme 
is R = log2 (NfM) bpcu. Thus, the SM scheme achieves 
an increase in spectral efficiency of logj Nt bits over single- 
antenna systems with a marginal increase in the complexity 
since it still needs only one RF chain at the transmitter 
However, in order to achieve high throughputs either Nt or 
M, or both need to be increased which renders this scheme 
suitable for low and moderately high spectral efficiencies. 

Recently, three specially tailored Sphere Decoding (SD) 
detectors were proposed for SM systems in JS] which were 
termed as Receiver-centric SD (Rx-SD) 14), Transmitter- 
centric SD (Tx-SD), and Combined-SD (C-SD) detectors. It 
was shown in ||3] that the Rx-SD detector is suitable for SM 
systems with large Nr, the number of receive antennas, and 
C-SD detector is suitable for systems operating at relatively 
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high spectral efficiencies i.e., with large Nt or M, or both. But, 
the applicability of the existing Tx-SD detector and hence, the 
C-SD detector is limited to systems with Nt < Nr, which is 
due to the zero diagonal entries in the R matrix of the QR 
decomposition associated with the underdetermined channel 
matrix. This problem is inherent to any MIMO system with 
underdetermined channel matrix that employs SD detector at 
the receiver. In ||5|, this problem was addressed by transform- 
ing the underdetermined channel matrix into a full-column 
rank matrix and applying standard SD on it. But, this SD 
which was termed as A-Generalized SD (A-GSD) detector is 
asymptotically Maximum Likelihood (ML) optimal with SNR 
for non-constant modulus signal sets such as 7\/-QAM. Thus, 
at low Signal-to-Noise Ratios (SNR) the performance of the 
A-GSD detector cannot be expected to be near-ML. The GSD 
of 161 is also not ML-optimal for non-constant modulus signal 
sets under all SNR conditions like the A-GSD detector, and 
also the complexity offered by the GSD of ||6l is significantly 
higher than that of the A-GSD fSl. In this paper, we do not 
take these approaches, instead, we show that the existing Tx- 
SD detector of [41 is not limited to SM systems with Nt < Nr 
but is applicable to systems with Nr < Nt < 2Nr — 1 as 
well, we refer to this detector as the Extended Tx-SD (E-Tx- 
SD) detector, and further, we propose an ML-optimal detection 
scheme termed as Generalized Tx-SD (G-Tx-SD) detector for 
SM systems with arbitrary Nt- 



II. System model 

We consider a MIMO system having Nt transmit as well as 
Nr receive antennas and a quasi-static, frequency-flat fading 
channel, yielding; 

y = Hx + n, (1) 

where x S C^*^^ is the transmitted vector, y G C^''^^ is 
the received vector, H G C^''^^' is the channel matrix, and 
n e ([^NrXi jjjg noise vector The entries of the channel 
matrix and the noise vector are from circularly symmetric 
complex-valued Gaussian distributions C7V(0,1) and CAf (0, 
a ), respectively, where ^ is the noise variance per dimen- 
sion. 



(2) 



A. Spatial Modulation 

In SM scheme HI, we have 

l-l Nt-l 

where s is a complex symbol from the signal set S with 
5*1 = M. Throughout this paper we assume 5* to be a lattice 
constellation such as QAM. Thus, for an SM system, eq.ldli 
becomes 

y = Hx;^s + n, (3) 

where I E L = and the subscript s captures the 

dependence of x on the signal set S. Assuming perfect channel 
state information and ML decoding at the receiver, we have 
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where yi and hi,i are the i* and the (/, i)* entry of the 
received vector y and the channel matrix H, respectively. 
The computational complexity in terms of number of real 
multiplications involved in computing ML solution of eq.(|5]l 
is given by 

(6) 
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since, 8 real multiplications are required in computing \yi — 
hi^is\'^ for any legitimate {s,l,i). If Nt > Nr, then from the 
well known division algorithm we have unique non-negative 
integers q and r (0 < r < Nr) such that Nt = qNr + r. Thus, 
from eq.® we can write 
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B. Review of Tx-SD detector ^ 

The complex- valued system in eq.([3j can be expressed in 
terms of real variables as 
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where, 5R(-) and represent the real and imaginary parts, 

y and n e M^jv.xi^ g ^ ^2N.x2Nt ^ and ^^^^ g Rsw^xi^ 

For Nt < Nr, the Tx-SD detector E) is given by 
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where, Qr — {{l,s) \ I G L, s e S, and ||y — Hx;^ 
with R as the initial search radius of the sphere decoder 
From ||4l and ||7|, we have R ~ aNra^, where a 
is a parameter chosen to maximize the probability 
of detection. By expressing H in terms of its QR 
decomposition we have ||y — Hx( ,j||2 < R^ equivalent 



to 



|z — Rx;.s||2 < Rq where, R is an upper triangular 



matrix given by 
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Fig. 1. Pictorial representation of tlie two level SD tree associated with 
each of the elements of x; The dependence of Xi on Xi_(-jVj for each i is 
indicated by grouping their branches in the same block. 



Q — [Qi {2N,.x2Nt) Q2 {2NrX2N^-2Nt)] SUCh that H — QR, 

and = i?2 _ ||Q|^y||2. Thus, we have 

= {(l.s) \leL,seS, and ||z - RiX/,,.||2 < R^q} ■ 

(10) 

Since Ri is upper triangular and x/ ^ has only two non-zero 
elements (from eq.® and eq.®) the elements of Qji are given 
by 

-Rq + Zi Rq+Zi 



< Xi < 
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' 2,2 

for iVt + 1 < i < 2Nt and 

— Rq + Zi — fi^i+Nt^i+Nt 



< X, < 



Rq + Zi — ri^i+N^Xi+Nt 



(12) 

for 1 < i < Nt, where, Xi and are the i* entries of x/ 
and z respectively, and fij is the entry of the upper 

triangular matrix Ri. From eqs. ( fTTT i and ( fT2b we have the SM 
specific SD tree shown in the Fig.® where N indicates the 
number of signal points in the real and imaginary dimensions 
of the constellation. For example, in a 64-QAM constellation, 
AI = 64 and iV = 8. 

III. Proposed Detector for SM systems with 

ARBITRARY Nt 

The Tx-SD detector discussed in the previous section was 
proposed for SM systems with Nt < Nr. In this section, we 
show that the existing Tx-SD detector is applicable to systems 
with Nr < Nt < 2Nr — 1 as well, and further extend it to 
systems with arbitrary Nt- 

A. Tx-SD detector in SM systems with Nr < Nt < 2Nr - 1 

Proposition 1 : The Tx-SD detector H originally pro- 
posed for SM systems with Nt < Nr is applicable to systems 
with Nr < Nt < 2Nr - 1 as well. 

Proof: Consider an SM system with 
Nt > Nr- The QR decomposition of H is 



given 
R = 



by QR where Q 
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where Ri is 



an upper triangular matrix. Recall that only two elements are 



non-zero in x/ and are apart by A^t — 1 zero elements as 
shown below. 
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Considering / = Nt and some s G 5 we have, 
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and the last element of p = RxjVt.s is given by P2Nr- = 
'R{2Nr, ■■)xNt,s- It is easy to see that R.{2Nr, :) has first 2Nr- 
1 elements zero since Ri is upper triangular, and if the number 
of non-zero elements 2Nt — 2Nr + 1 in R(2A',., :) is less than 
or equal to Nt = {Nt — 1) + 1, the number of zeros between 
5R(s) and 3(s) plus one non-zero element 3(s) in i<-Nt,s, we 
see that p2Nr depends only on 5J(s). Since this is true for all 
{Pi}'i=2 ™y legitimate pair {I, s), we have the condition 
2Nt - 2Nr + 1 < Nt = Nt < 2Nr - 1 for independent 
detection of imaginary components in x; s. This results in the 
following intervals analogous to those in eq.lfTTI). 
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'''2N,^-i+l,2Nt-i+l r2Nr.-i+l,2Nt-i+l 

for I < i < Nt- Proceeding in the Unes similar to that of 
im for real components in x; s we get the intervals of eq.(fT6]l 
given in the next page. From eq.(fT6b we observe that x. With 
Nt = Nr, it can be checked that the intervals of eq.lfTSl) and 
eq.lfTSb reduce to those of eq. lfTTl i and eq. (fT2l) . respectively. ■ 
We refer to this Tx-SD detector as the Extended Tx-SD (E- 
Tx-SD) detector in the rest of the paper 

B. Tx-SD detector for SM systems with arbitrary Nt 

We have shown in the previous subsection that the Tx-SD 
detector of Q can be used in SM systems with Nr < Nt < 
2Nr — 1. In this subsection we propose a SD detection scheme 
for arbitrary Nt by partitioning the antenna search space into 
disjoint subsets each of size 2Nr — 1 and running E-Tx-SD 
decoders sequentially. 

Consider an SM system with Nt > 2Nr - 1 ^ N'. From 
the division algorithm we have unique non-negative integers 
q and r such that Nt = qN' + r, where < r < iV'. For the 
ease of presentation we assume r = and hence Nt = qN' 
for now, and later generalize our results for non-zero r. Let L, 
the set of antenna indices, be partitioned into disjoint subsets 

Lk = {«}*^=('fc_i)jv'+i so that L = ULi^fc and \Lk\ = N' 
for all 1 < fc < g. 

Let J{l,s) represent the ML metric X^il^i IVi ~ hi,is\'^ of 
eq.(|5]i. Then, the ML solution in terms of J{1, s) is (l, s)ml = 
argminjgi J(Z, s) and we have 
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\Je{k) , (19) 

where 8(fc) = minjgi^. {min^gg J(/, s)}. In the above, 
(a) follows from the assumption that the antenna index and 
the transmitted symbol are encoded by independent sets of 
bits, and (b) follows directly from the fact that L partitions 
into Lfc's. It is straightforward that each of the 0(fc)'s can 
be obtained by running E-Tx-SD detector discussed in the 
previous subsection. The sphere radius Rk for all 1 < fc < g 
is taken as i? = aNrCr"^- Thus, we have 

Q{k) = {{l,s) I ?eifc,se5, and ||zfe - R^x^.^H^ < i?U , 

(20) 

where, QfcRfc = H(:, [Lfe]), = R^ and Zfc = Q^y 

for all 1 < fc < (J. Equations (flSi l and ( fT6] l directly give the 
elements of 0(fc) for all 1 < fc < 

Now, for a system with Nt = qN' + r and r 7^ 0, it is 
straightforward that there is an additional set 0(5 + 1) in 
eq.(fT9]l which is same as eq.(l20li with i^+i ^ L \ [Jl^i Lk- 
Thus, the E-Tx-SD based solution for systems with arbitrary 
Nt is given by 
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This solution is referred to as the Generalized Tx-SD (G-Tx- 
SD) detector. We note here that the partitions Lk considered 
here are of size N' ~ 2Nr — 1, it is straightforward that the 
detector in eq. (l2Tl l can be run over partitions of any size, for 
example N' = Nr as well. 

From eq.© it is clear that the ML complexity scales linearly 
with the number of transmit antennas. As the proposed G- 
Tx-SD detector runs over partitioned channel blocks, it is 
obvious that the complexity of the proposed detector can be 
expected to be lesser than that of the ML detector, since the 
complexity of the individual Tx-SD/E-Tx-SD detector is much 
lesser than that of the ML detection |4|. We note here that the 
proposed E-Tx-SD and G-Tx-SD detectors are not limited to 
SM scheme alone, but are applicable to any system with the 
ICl-free property, for example, the Space-Time Shift Keying 
ID. The C-SD detector proposed in [4] uses both the Tx-SD 
and the Rx-SD detectors in order to reduce the overall search 
complexity. It is straightforward that the proposed detectors 
can be further extended by incorporating Rx-SD detector for 
reduction in the receive search space complexity as well. 

IV. Simulation results 

Consider an SM system with Nr — 4: and Nt = 7, 
employing 16- and 64-QAM signal sets. Let Nml ~ MNt 
denote the number of search points in the ML detection, and 
¥,[Ne-Tx-sd] denote the expected number of points that lie 
inside the hypersphere of radius R under E-Tx-SD detection, 
in this paper, we use the metric K[Ne-Tx-sd]/Nml to 
measure the reduction in search space with respect to the 
ML detection. From Fig.®, it is clear that the Symbol Error 
Rate (SER) performance of the ML and the E-Tx-SD detectors 
overlap for both the modulation schemes considered, and from 
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Fig. 2. SER curves of ML and E-Tx-SD detectors in an SM system with 
A^r = 4:, Nt = 7, employing 16- and 64-QAM signal sets. 

Fig.© we see that the expected number of nodes that lie 
inside the sphere reduces significantly with increase in SNR. 
Specifically, at an SNR of 20 dB, we see a reduction of about 
90% and 70% for 16- and 64-QAM signal sets, respectively. 
Fig.® shows the SER performance and the reduction in 
complexity due to the G-Tx-SD detector in an SM system 
with Nr = A, Nt ^ 8, and 16-QAM signal set. We observe 
from Fig. ©(a) that the SER curves of the proposed and the 
ML detector overlap, and from Fig.(|4|(b) that there is about 
85% reduction in complexity with respect to ML detection at 
an SNR of about 18 dB. The partitions considered here are of 
size \Li \ = 1 = 4. The reduction in complexity involved in 
precomputations such as QR decomposition, etc., by using the 
E-Tx-SD detector with different |Li|'s instead of the Tx-SD 
detector with equal \Li\'s is left for our future study. 

V. Conclusion 

We have shown that the existing Tx-SD detector is not 
limited to SM systems with Nt < N,- but is applicable to 
systems with A^,. < Nt < 2Nr — 1 as well. Further, we 
have proposed a generalized Tx-SD detector for SM systems 
with arbitrary Nt by exploiting the ICI-free nature of the 
system, and have shown with our simulation results that the 
proposed detectors give ML performance with significantly 
reduced complexity. 
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